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We confront the most common version of the CP-conserving 2HDM with LHC data, taking into 
Q 



account all previously available experimental data. We also discuss the scenario where the 125 GeV 
Higgs discovered at the LHC is the lightest neutral scalar of a particular CP-violating 2HDM. In this 

Oh scenario we focus on what data can already tell us about the amount of mixing between CP-even 

D and CP-odd states. 



^ I. INTRODUCTION 

(N 

C\| After the discovery [TJ[2] of a Higgs boson at CERN's Large Hadron Collider (LHC), it is now time to confront 

possible extensions of the Standard Model (SM) with LHC data. One of the simplest ways to extend the scalar 
■^j- sector of the SM is to add one more complex doublet to the model. The resulting two Higgs doublet models 

(2HDMs) can provide additional CP- violation coming from the scalar sector and can easily originate dark matter 
candidates. More evolved models with additional field content have a 2HDM like scalar sector - as in the case, 
for instance, of the Minimal Supersymmetric Standard Model. 2HDMs have a richer particle spectrum with one 
£> charged and three neutral scalars. All neutral scalars could in principle be the scalar discovered at the LHC [3]- 
[5]. However, data gathered so far allow us to rule out the 125 GeV particle as being a pure pseudoscalar state. 
Hence, we are left with either one of the CP-even states of a CP-conserving model or a neutral state from a 
CP-violating model [BJ. In the next section we will present the 2HDM models and in the following sections we 
discuss in turn how well 2HDMs can accommodate LHC data for the CP-conserving and for the CP-violating 
case. 

II. THE MODELS 

In this section we present the two versions of the 2HDM to be discussed, one CP-conserving and the other 
explicitly CP-violating. The most general 2HDMs' Yukawa Lagrangian gives rise to scalar exchange flavour 
changing neutral currents (FCNCs) which are strongly constrained by experiment. A simple and natural way 
to avoid those dangerous FCNCs is to impose a Z^ symmetry on the scalar doublets, $i — > $i, $2 — > —$2, and 
a corresponding symmetry to the quark fields. This leads to the well known four Yukawa model types I, II, Y 
(III, Flipped) and X (IV, Lepton Specific) [TJ. [B] . The different Yukawa types are built such that only </> 2 couples 
to all fermions (type I), or <f>2 couples to up- type quarks and <f>\ couples to down- type quarks and leptons (type 
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II), or fa couples to up- type quarks and to leptons and fa couples to down- type quarks (type Y) or finally fa 
couples to all quarks and fa couples to leptons (type X). 

The scalar potential in a softly broken Z 2 symmetric 2HDM can be written as 

V(*i, $ 2 ) =77-1?$!$! + m|*£*2 + (m 2 2 $j$ 2 + h.c.) + ~Ai(*I*i) 2 + ^A 2 (^$ 2 ) 2 

+A 3 ($ t 1 $ 1 )($|$ 2 ) + A 4 ($ t 1 $ 2 )($ t 2 $ 1 ) + ~A 5 [($1$ 2 ) 2 + h.c] , 

where $j, i = 1, 2 are complex SU(2) doublets. All parameters except for m\ 2 and A5 are real as a consequence 
of the hermiticity of the potential. The phases of m 2 2 and A5 together with the ones in the vacuum expectation 
values (VEVs) will determine the CP nature of the model (see [5] for a review) . A CP-violating model has three 
spin neutral states usually denoted by h\ , h 2 and h 3 while in the CP-conserving case the three CP-eigenstates 
are usually denoted by h and H (CP-even) and A (CP-odd). As shown in [TOj, [TT] , once either a CP-conserving 
or a CP-violating vacuum configuration is chosen, all charge breaking stationary points are saddle points with 
higher energy. Hence, the 2HDM is stable at tree-level against charge breaking and once a non-charge breaking 
vacuum is chosen both models have two charged Higgs bosons that complete the 2HDM particle spectrum. 

We will focus on two specific realisations of 2HDMs. One is the usual 7-parameter CP-conserving potential 
where m 2 2 , A5 and the VEVs are all real. In this model we choose as free parameters, the four masses, the 
rotation angle in the CP-even sector, a, the ratio of the vacuum expectation values, tan/3 = v 2 /v\, and the soft 
breaking parameter redefined as M 2 — m 2 2 /(sin/3 cos/3). The other is an explicit CP-violating model [12l414j . 
In the CP-violating version m\ 2 and A5 are complex and the fields' VEVs are real. Existence of a stationary 
point requires Im(m 2 2 ) = v\v 2 Im(\§). Because the VEVs are real in both models, a common definition for 
the rotation angle in the charged sector tan/3 = v 2 jv\ can be used. In this model, besides mjy± and tan/3, we 
take as free parameters the masses of the lightest scalars m^i and rtihi the three rotation angles in the neutral 
sector «!, a 2 and a.3 and Re(m 2 2 ). Details of the model can be found in [T2TrH] . 

A. Constraints on the models 

We have imposed the following theoretical bounds on both models: we require that the potential is bounded 
from below |15j and we impose unitarity limits on the quartic Higgs couplings 16J. We have also taken into 
account the precision electroweak constraints [TTJ HH]. Furthermore, we discuss the scenarios in which two 
normal minima coexist in the CP-conserving potential. In those scenarios it could happen that we are living in 
a local minimum and tunnelling to a global deeper minimum could occur. This deeper minimum would have a 
different v 2 and therefore different masses for the elementary particles. We call this a panic vacuum [TTI [T9"l |2"U] . 
The experimental bounds pertaining to the charged sector apply to both models since the charged Higgs Yukawa 
couplings have exactly the same form. Except for the LEP bounds they all constrain regions of the (tan /3, m,H± ) 
plane. The LEP experiments have set a lower limit on the mass of the charged Higgs boson of 80 GeV at 95% 
C.L., assuming BR(H+ -> t+v)+BR{H+ -> cs)+BR(H+ -> AW+) = 1 |H] with the process e+e~ -> H+H~. 
The bound on the mass is 94 GeV if BR(H + — >• t + v) = 1 [5T]. These are model independent bounds as long as 
the above mentioned sum of BRs is one. Values of tan/3 smaller than O(l) together with a charged Higgs with a 
mass below 0(100 GeV) are both disallowed by the constraints [22 coming from R^, from B q B q mixing and from 
B — > AT S 7 [33] for all models. Furthermore, data from B — > X s j [33] imposes a lower limit of m H ± > 360 GeV, 
but only for models type II and type Y. We have also considered the most recent bounds from the ATLAS [24! 
and CMS [25] collaborations on the (tan/3,m#±) plane coming from pp — > tt — > bbW T H ± (— > tv). Finally, we 
have considered the LEP bounds on the neutral Higgs 26 . 

Combining it all, our scans will be performed taking m H ± > 90 GeV and tan /3 > 1 for type I while m H ± > 360 
GeV for type II. For the CP-conserving case we also take vha > 90 GeV, rrih = 125 GeV and tuh > m^ while a 
is free to vary in its allowed range although subject to the above constraints. For the CP-violating model we take 
rri/ji = 125 GeV and m/ l2 > mhi while on [i — 1, 2, 3) are again free to vary in their allowed ranges, subject to 
the above constraints. Finally, we note that we do not include the constraint that would result from the anomaly 
observed by the BaBar collaboration in the rates R(D) and R(D*), with R(D) = (B — >• Dt~V t )/(B — > Dl~Vi), 
which deviates by 3.4 a (when D and D* final states are combined) from the SM prediction [27 . It is clear that 
the deviation cannot be explained by a charged Higgs boson from either of the four FCNC conserving Yukawa 
types. Therefore, if this observation is independently confirmed by the BELLE collaboration not only the SM 
will be excluded at 3.4 a but all 2HDMs Yukawa types will also be excluded with a very similar significance. 

We should point out that there will never be an exclusion of the 2HDM in favor of the SM. Indeed, the models 
discussed here include the SM as a specific (decoupling) region of parameter space. As a result, experiments 
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may exclude the SM in favor of the 2HDM; they may exclude both; or they may be compatible with SM, thus 
restricting the 2HDM extension into a particular (eventually, small) region of the parameter space. 



III. CP-CONSERVING MODEL 



We start with the CP-conserving model and randomly generate points in the parameter space such that 
m h = 125 GeV, 90 GeV < m A < 900 GeV, m h < m H < 900 GeV, 1 < tan/3 < 40, -(900) 2 GeV 2 < m\ 2 < 900 2 
GeV 2 and -tt/2 < a < tt/2. Further, for type I we take 90 GeV < m H ± < 900 GeV while for type II the 
allowed range is 360 GeV < m H ± < 900 GeV. In order to use the LHC data we define the quantities Rf as the 
ratio of the number of events predicted in the 2HDM to that obtained in the SM for a given final state /. 



Rf = 



OJPP -> fe) 2 HDM BR(fe -» /) 2 HDM 

a(pp -> h SM ) BR(/i SM -> /) 



(1) 



where h is the lightest CP-even Higgs (125 GeV), a is the Higgs production cross section, BR the branching 
ratio, and /ism is the SM Higgs boson. In our analysis, we include all Higgs production mechanisms, namely, 
gluon-gluon fusion using HIGLU at NLO [28], vector boson fusion (VBF) [29], Higgs production in association 
with either W, Z or tt [25], and bb fusion [315]. A different approach was recently followed in [3T] where they 
have looked for the best fit point in parameter space for the four Yukawa types of the same CP-conserving 
2HDM. 
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FIG. 1: Points in the (sina, tan/3) plane that passed all the constraints in type I (left) and in type II (right) at la in 
green (light grey) and 2<r in blue (dark grey). Also shown are the lines for the SM-like limit, that is sin(/3 — a) — 1 
(negative sina) and for the limit sin(/3 + a) = 1 (positive sina). 



The points generated have then to pass all the constraints previously described plus the combined [32] ATLAS 
and CMS strengths i? 77 = 1.66 ± 0.33, R ZZ = 0.93 ± 0.28 and R TT = 0.71 ± 0.42. It is important to point out 
that this talk was given in February, 2013, prior to the updates from Moriond. We have added a subsection 
regarding the Moriond updates below. In figure fl] we present the points in the (sina, tan/3) plane that pass all 
the constraints in type I (left) and in type II (right) at la in green and at 2a in blue. Also shown are the lines 
for the SM-like limit, that is sin(/3 — a) = 1 (negative sina) and for the limit sin(/3 + a) = 1 (positive sina). 

To better understand these results we approximate Rzz in the case where Higgs production is due exclusively 
to gluon-gluon fusion via the top quark loop, and the total Higgs width is well approximated by T(h — ¥ bb). 
This is a very good approximation for most of the parameter space and it yields 



R 



I,approx 

zz — 



cos a 
sin 2 /3 



sin 2 (a — /3) 



sin 2 /3 



= sin 2 (a — /3). 



(2) 



By setting R£ z prox = 1 one obtains the SM-like limit line shown in the picture for negative a. Applying the 
same simplified scenario we obtain for type II 



R 



1 1 ,approx 
ZZ 



(a-/3) 



1 



tan a tan j3 



(3) 



Toyama International Workshop on Higgs as a Probe of New Physics 2013, 13-16, February, 2013 



II,approx 

ZZ 



1 we obtain not one but two lines, the red lines shown in figure 1 The approximate 



ibo 



Now, when R 

scenarios lead us to the following conclusions. First, it is very hard in type I to have Rzz ab~ove 1. Further, the 
function sin (a — j3) is very sensitive to deviations from 1. This leads to a larger dispersion of points around 
the SM-like limit line. Regarding type II, not only do we get two lines instead of one but the function R z ^ pprox 
is now very insensitive to deviations from 1 and that is why the points are all in a narrow band close to each 
line. It is interesting to note that while the tension between Rzz and -R 77 excludes all points at la in type II, 
there is still a band in type I close to sin a = where Rzz is below 1 while i? 77 can be large. The la points in 
type I all have a charged Higgs mass below 130 GeV. 
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FIG. 2: In the left panel we show the points in the (sina, tan/3) plane that have passed all the constraints in type II 
with an exact Zi symmetry. In the right panel we show type II with the broken Zi symmetry where we have plotted 
the panic points in red. No points survive at la and the 2a points are shown in blue. 



Type I - Decoupling Limit 




Type II - Decoupling Limit 




FIG. 3: Points in the (sina, tan/3) plane that passed all the constraints in type I (left) and in type II (right) at 2<r in 
blue. We have taken all masses except rrih = 125 GeV to be above 600 GeV. 

In the left panel of figure [2] we present the points that have passed all the constraints in type II with an 
exact Z 2 symmetry {m\ 2 = 0). It was recently shown [33] that if rrih = 125 GeV, then tan/3 has to be below 
approximately 6 . This is confirmed in the plot. Apart from the limit on tan/3, the plot does not differ much 
from the softly broken Z 2 case except that there are fewer points and that they move slightly away from the 
red lines when compared with the soft breaking scenario. We note that in type I no points survive at la in the 
exact Z 2 symmetric case. In the right panel of figure [2] we show the panic points (in red) for the soft breaking 
type II model. Panic points are easily found by using the following discriminant [2D] 



D = (m?! - k 2 m\ 2 ) (tan/3 - k). 



(4) 
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where k — yAi/A^- li D < we live in a metastable state and our current (t>i, V2) solution is the panic vacuum. 
We found panic points distributed all over the 2a allowed region. We should note however that as shown in |20j . 
the panic points are further away from the SM-like limit than the non-panic points. 

In figure [3j we present similar plots to the ones shown in figure IT] the only difference being that all masses 
except for the lightest Higgs one are taken to be above 600 GeV. This plot is representative of the decoupling 
limit of the 2HDM - if all particles except h are too heavy to be detected at the LHC this is how the plot would 
look after the bounds on those masses are considered. 

We finish this section with a comment on the values of M 2 and the values of the remaining masses that 
are still allowed by data. We had already noted that the points with M 2 = are disallowed at 2a in type 
I. We have checked that also M 2 < is excluded at 2a in type I. In type II all values of M 2 are allowed at 
2a. Regarding the remaining masses, our conclusion is that the data on the 125 GeV Higgs does not provide 
meaningful constraints on the masses of the remaining scalars. 



A. Update after Moriond 2013 

In this section we present updated plots after Moriond 2013. We use the latest updates from from ATLAS [51] 
and from CMS [35] for the 7 and 8 TeV runs. Regarding the CMS results for h — > 77 we use the result from 
the Multivariate analysis. 
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FIG. 4: Points in the (sin a, tan/3) plane that passed all the constraints in type I using the ATLAS data (left) and using 
the CMS data (right) at la in green (light grey) and la in blue (dark grey). Also shown is the line for the SM-like limit, 
sin(/3 — a) — 1. 



In figure|4]we present the results for the type I model in the (sin a, tan/3) plane using the ATLAS data (left) 
and the CMS data (right). The differences between the two plots are easy to understand. The ATLAS data 
forces Rzz to be large but as we have previously shown (equation vM Rzz can never be above one in type I. 
Consequently, no points survive at la. With the CMS data plenty of la points survive because all Rvv are 
below one. As expected, the la region is slightly away from the SM-like limit because, on one hand, the central 
values of Rvv are below one, and, on the other hand, we already saw that Rvv are very sensitive quantities in 
type I. 

In figure [5] we now show similar plots for the type II model. In type II Rvv is not a sensitive quantity. 
Therefore both the la and the 2a points are very close to the two limiting lines (in red) . The difference between 
the two plots is that the points tend to concentrate below the limiting lines for ATLAS and above the same 
lines for CMS. This is of course a consequence of the central values of the ATLAS Rvv being above the SM 
while the CMS ones are below the SM expectation. Finally we present in figure [6] plots similar to the ones in 
figure p\ but now in the (cos(/3 — a), tan/3) plane. Again we see the same trend with the allowed points on 
opposite sides of the limiting lines, cos(/3 — a) = and cos(/3 + a) = 0. 
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FIG. 5: Points in the (sin a, tan/3) plane that passed all the constraints in type II using the ATLAS data (left) and using 
the CMS data (right) at la in green (light grey) and 2a in blue (dark grey). Also shown are the lines for the SM-like 
limit, that is sin(/3 — a) = 1 (negative sin a) and for the limit sin(/3 + a) — 1 (positive sin a). 
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FIG. 6: Points in the (cos(/3 — a), tan/3) plane that passed all the constraints in type II using the ATLAS data (left) 
and using the CMS data (right) at la in green (light grey) and 2a in blue (dark grey). Also shown are the lines for the 
SM-like limit, that is cos(/3 — a) = and for the limit cos(/3 + a) — 0. 



IV. CP- VIOLATING MODEL 



In this section we discuss the CP- violating model 2HDM presented in section [TTJ The lightest neutral particle 
of the model h\ is considered to be the one discovered at the LHC. This particle can be seen as being composed 
by a CP-even and a CP-odd component. The amount of mixing is controlled by S2 = sina2 (one of the three 
rotation angles in the neutral sector - see [6] for details). When S2 — 0, hi is a pure CP-even state and when 
,s 2 = 1, hi is a pure CP-odd state. Our goal is to constrain the amount of mixing using LHC data. In the left 
panel of figure [n we show Rzz versus _R 77 for type I and for three regions of the parameter s 2 : \s%\ < 0.1 in 
green (light grey), 0.45 < |s 2 | < 0.55 in blue (black), and |s 2 | > 0.83 in red (dark grey). In the right panel 
we present a similar plot for type II. We also show in both plots the combined ATLAS and CMS results at 
lcr (yellow line) and at 2a in blue. We conclude that the red points are excluded at 2a when the combined 
results [32] from ATLAS and CMS are considered. 

The results of the Moriond updates are shown as 2a solid (CMS) and dashed (ATLAS) lines in figure [7| We 
see that | s 2 1 > 0.83 is still excluded at 2a by the ATLAS data, but that it is allowed by the new CMS data. 
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FIG. 7: In the left panel we show Rzz versus R 11 for type I and for three regions of the parameter S2' |s2| < 0.1 in 
green (light grey), 0.45 < \sa\ < 0.55 in blue (black) and | «2 1 > 0.83 in red (dark grey). In the right panel we present a 
similar plot for type II. We show in both plots the combined ATLAS and CMS results at la (yellow line) and at 2a in 
blue and the 2a updated lines after Moriond 2013 (solid black - CMS and dashed black - ATLAS). 

V. CONCLUSIONS 

We have considered the most common CP-conserving 2HDMs and a CP-violating 2HDM in light of the recent 
data from the LHC. For the CP conserving case, there are 4 different models, corresponding to the transformation 
properties of the quarks, leptons and Higgs bosons under a Z 2 symmetry. We include constraints from vacuum 
stability, unitarity, precision electroweak constraints, LEP experimental bounds, bounds from b decays as well 
as LHC bounds on top decays. The full parameter space is scanned, subject to these constraints, and allowed 
points are plotted in the (sin a, tan (3) plane and in the (sin a, cos(/3 — a)) plane. This talk was presented before 
the Moriond 2013 meeting, at which time the Standard Model branching ratios disagreed by more than one 
sigma with the combined data from ATLAS and CMS, and thus the plots show no points within one sigma 
of the Standard Model curves. However, we have updated the plots since Moriond, and find that in type I 
CMS data does allow points within one sigma of the Standard Model. In type II both ATLAS and CMS now 
allow points within one sigma of the SM. In the next run, these experimental bounds will tighten considerably, 
providing a strong test of the models. 

In the CP violating case, the key parameter is S2, which controls the amount by which the observed 125 GeV 
state is CP-even vs. CP-odd. We have shown that some regions of relatively large s 2 can be excluded, but a 
wide range of values remains. Over the next few years, the rectangles shown in figure [7] will shrink considerably, 
providing a much stronger bound on the amount of the scalar's CP-odd component. 
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